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Abstract 

We report on magneto-transport measurement in disordered graphene under pulsed magnetic 
field of up to 57T. For large electron or hole doping, the system displays the expected anomalous 
Integer Quantum Hall Effect (IQHE) specific to graphene up to filling factor v = 2. In the close 
vicinity of the charge neutrality point, the system breaks up into co-existing puddles of holes and 
electrons, leading to a vanishing Hall and finite longitudinal resistance with no hint of divergence 
at very high magnetic field. Large resistance fluctuations are observed near the Dirac point. They 
are interpreted as the the natural consequence of the presence of electron and hole puddles. The 
magnetic field at which the amplitude of the fluctuations are the largest is directly linked to the 
mean size of the puddles. 



PACS numbers: 81.05.Uw , 71.20.-b 



A few years ago, a new form of IQHE has been discovered in graphene for which the 
Hall conductance is half-integer quantized in units of 4.e 2 /h l|. The factor 4 stands for 
Landau level (L.L.) degeneracy (two-fold spin and two-fold nodal degeneracy) so that 
IQHE occurs at filling factor v = 4 x (|z| + 1/2), where i — 0, ±1, ±2... is an integer used 
for L.L. indexation. In very high magnetic field and close to the Dirac point, spin and 
nodal degeneracies are expected to be lifted. Indeed, pioneering experimental works 
showed additional Hall resistance plateaus for filling factor v — 0, v — ±1, v — ±4 as well 
as ±3. While Zeeman coupling lifts spin degeneracy, the nodal degeneracy lifting implies 
spontaneous svmmetry breaking effects 4j, |5| which fundamental nature is still under 
investigation |6h12|. It is now established that disorder significantly affects the electronic 
properties of graphene and thus cannot be ignored, especially at low charge carrier density 
where it brings about the formation of electron-hole puddles Jl3|. Even if the field-induced 
insulator transition is well documented for clean graphene [14J, |l5|, little is known about 



high-field transport in disordered samples where the pres ence of electron-hole puddle is 
expected to dominate over a large magnetic field range 



In this communication, we report on IQHE in a disordered system for various charge 
carrier densities from hole to electron doping. At low field, our experimental data are 
consistent with existing theories for disordered graphene. Close the the Dirac point, the 
magneto-resistance is quasi-flat for a broad magnetic field range without sign of divergence 
in contrast with the clean case. In addition, at large magnetic field, we show that both 
longitudinal and Hall magneto-resistance display large fluctuations at the Dirac point, whith 
a maximum amplitude for B m 25T. We show that these fluctuations result from the 
presence of electron-hole puddles with an estimated mean diameter of 5nm. 

The sample has been designed in the standard Hall bar geometry using electron-beam 
lithography, allowing simultaneous measurements of both Hall (R X y) and longitudinal {R xx ) 
resistances. Variable electrostatic doping is achieved upon applying a voltage through 
300w7i of S1O2 dielectric used as a back-gate. Figure [TJ shows the device resistance as the 
gate voltage V g is ramped from — 70V to +70V. The charge neutrality point is identified 
at the resistance maximum Vqnp — +52V. The large FWHM of the curve as well as 
the large shift of Vcnp from V g = 0V suggests a rather disordered sample. Actually, 
non-intentional hole doping certainly arises from defects and charged chemical species 
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FIG. 1: Longitudinal resistance as a function of gate voltage V g at T = 1.5K and B = OT, the 
charge neutrality point is located at V g = +52V. The arrows indicate specific values of gate voltage 
at which magneto-resistance data will be presented (see figure [3]). Inset : Longitudinal resistance 
R xx and Hall resistance R xy as a function of magnetic field at T = 1.5K and V g = OV. 

adsorbed at the surface of the sample or at the S1O2 interface during its fabrication or 
exposure to air. The gating efficiency of the device is estimated to /3 = 7.85 x 10 10 cm~ 2 /V 
from low field Hall measurements. This result is compared to the accepted plane capacitor 
model for graphene = (e e r )/(e.cf) where e$e r is the dielectric permittivity of Si02 and 
d = 275nm its thickness, in good agreement with the nominal value. A constant current of 
i = 0.5fiA is injected through the device at 1.5K while a strong pulsed magnetic field, of 
total duration 300ms and rising up to 57T, is applied perpendicular to the sample. Inset 
of figure [T] shows IQHE at V g = 0V. The Hall resistance displays well defined quantized 
plateaus at R xy = + 1/2) -1 for i — 0, 1, 2... typical for graphene. The mean mobility is 
estimated to /i = 1300cm 2 /V!s in consistency with the first appearance of SdH oscillations 
at minimum magnetic field m IT (see inset of figured]). 

For B < B>, in the classical regime, the experimental data can be approached using a 
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FIG. 2: (a) and (b) Theoretical adjustment of R xy (B) and R xx (V g ) at small field respectively, at 
selected back gate voltages or magnetic field using the two-fluid model. R xx (V g ) data have been 
successively upshifted of AR = 5A;$7 for clarity. Note that at B = B] = 7T, the two-fluid model is 
no longer accurate (c) and (d) Extracted coefficient a and mobility \i as a function of V g : notice 
the coexistence of electrons and hole-like particles for a broad range of gate voltage AV g ~ 23V, as 
well as a mobility dip near CNP. Dotted lines are guide for the eye. 



standard two-fluid model 17|, |l8j. In this picture, electron and hole-like particles coexist and 
both contribute to transport. For simplicity and within a good approximation, we assume 
the same mobility [i for both electrons and holes. Their density will be denoted n e and 
respectively. In the following, we use the reduced parameter a = (rih — n e )/(n e + nh) so that 
R XX (B) = R xx (0) x and R xy = ^a.^.B.R xx (B) with R xx (0) = 7 x [(n e + n^.e.fi]- 1 



-(a.fi.B) 

and 7 = 3 is the form factor |19| . A theoretical adjustment with a and \x as free parameters 
is performed to match the full experimental dataset. The agreement is convincing as shown 
in figure [2]-a and |2]-b. We emphasis on the fact that R xx (0) is not a free parameter and is 
determined from the R xx (V g ) curve at zero magnetic field (see figure d]). Although some 
refinements (e.g. different mobilities for electrons and holes particles) would certainly 
improve the fitting quality, both the longitudinal and Hall resistance can be adjusted with 
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this minimal two-fluid model which captures well the underlying physics of disordered 
graphene. This fitting procedure allows us to extract the profiles of ot(V g ) and n(V g ) (see 
figures |2]-c and[2]-d). On crossing the charge neutrality point, the hole to electron density 
ratio changes smoothly within a gate voltage interval AV g ~ 20V. In the presence of 
disorder, mainly modeled by charged impurities, the coexistence of hole and electron-like 
particles is expected within a window AV g = where n* is an effective charge carrier 

(5x10 15 /mT 



density defined by n* 



I.P.VCNP 



20|. We find AV g ~ 23V in good agreement with 



experimental results. The predicted minimum of conductivity at CNP and zero field is 
Gcnp = hj x x 5*10^ = 61/iS 1 which is favorably compared with experiment (i.e. 
G £ Jnp — 60.6/iS 1 ). In first approximation, graphene's resistivity was experimentally found 
to be inversely proportional to the density of charge carriers, so that their mobility is 
almost independent of n oc V g . This particularly holds for low mobility samples, far from 
the Dirac point, where conduction is mainly dominated by long-range scatterers such as 
charged impurities. On the other hand, little is known about the carrier's mobility near 
the CNP since, when calculated using simple Drude model, fi = a/e.n diverges for n — > 0. 
Using the two-fluid model fitting procedure, the extracted curve n(V g ) in figure [2]-d shows 



21| compared 



a clear minimum of the mobility in the vicinity of the CNP. Do and Dollfus 
the effects of long-range and short-range scattering centers on mobility. They demonstrated 
that, even if the mobility is peaked at V g — Vcnp i n the presence of charged impurities, it 
is strongly reduced when vacancies are included, leading to a dip-like curve consistent with 
our observations. 



We now turn our attention to high magnetic field transport properties. Figure [3]-a and 
[3]-b show the longitudinal and Hall resistances in the vicinity of the charge neutrality point. 
Focusing close to Vcnp ~ 52V, the longitudinal magneto-resistance is first positive and 
reaches a maximum for B « B' = 7T. This resistance maximum is attributed to the forma- 
tion of overlapping LLs and naturally sets the crossover from the classical to the quantum 
regime. The resistance then decreases as a consequence of increasing LLs energy interval 
with magnetic field till the i = LL only remains populated. From B f» 20T, a merely flat 
magneto-resistance is observed up to the highest experimentally available magnetic field. 
In graphene, the energy of the n = LL is independent of magnetic field. Therefore, as 
long as long as the Fermi energy remains at a constant value within the n = Landau 
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FIG. 3: Longitudinal resistance (a) and Hall resistance (b) as a function of magnetic field for selected 
gate voltages (see arrows in figure [1]) in the vicinity of the charge neutrality point. Note that in 
(a), the curves have been successively down-shifted of — 5k$l with respect to curve at V g = 30V for 
clarity. Inset shows R XX (B = 40T). 

Level, a fixed back-scattering rate between edge states through extended states accounts 
for the quasi-flat magneto-resistance. The scattering rate is maximum when V g = Vcnp, 
that is when the Fermi energy is aligned with the center of the n = L.L. (see inset of 
figure [3]). As a matter of fact, instead of a purely flat magneto- resistance, a weak increase 
is observed at very high magnetic field. The closer V g is to Vcnp, the less pronounced is the 
high field magneto-resistance increase. This effect is explained considering a Fermi energy 
shift towards the center of the n = LL as its degeneracy increases with the field. The 
energy shift of Ep is all the more important as the initial density of states is weak. We 
note that, contrary to clean graphene devices [22], no divergence of the magneto-resistance 
is observed up to B = 57T for this particular disordered device. This figure, together with 
the above-mentioned R xx (B,V g ) behavior, rule out any alternative explanations implying 
spontaneous symmetry breaking and degeneracy lifting. However, whether a divergent 
behavior will be detected or not for even higher magnetic field remains an open question. 
In the close vicinity of the CNP, the Hall resistance displays no more well defined quantized 
plateaus (e.g. compared with R xy (B) curves at V g = 30V showing well developed Hall 
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plateau at R xy = 12.9&0), but remains close to zero as a result of electron and hole balance. 
We note that, for some fixed values of the gate voltage, the Hall resistance experiences a 
sign inversion indicating a change of the majority carrier type. This effect is related to a 
field-enhanced coexisting energy zone of electron and hole particles and will be discussed 
elsewhere. 




FIG. 4: (a) and (b) : Longitudinal and Hall resistance fluctuations for various gate voltages. A 
smooth background function has been subtracted to the original data. Data have been shifted 
for clarity, (c) Standard deviation of the curves shown in a) (black) and b) (red), d) Maximum 
peak-to-peak amplitude of the fluctuations. 

In the vicinity of the Dirac point, the Hall and longitudinal magneto-resistances are 
distorded with large fluctuations. Such fluctuating features are qualitatively analyzed in 
the conductance data, obtained by resistance tensor inversion. Figures 0]-a) and HJ-b) high- 
light on such conductance fluctuations, where a smooth background function G(B) has 
been subtracted to the original conductance data. The standard deviation as a function of 
gate voltage has been displayed in figure H]-c). The amplitude of the fluctuations are much 
larger in the vicinity of the CNP (figure HJ-c) and in the intermediate magnetic field range 
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(10T < B < 35T). Similar fluctuations have already been reported in various DC high 
magnetic-field experimental works j^, [ll| as a function of gate voltage, but they are not 



systematically observed 



23fl depending on sample quality, magnetic field range and tem- 



perature. In disordered materials, magneto-resistance fluctuations are usually understood 
in terms of magnetic-field dependent phase shifts of the electronic wave-functions and are 
well reproduced using the UCF theory 
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25] . However, their occurrence in the vicinity 
of the Dirac point only and their large amplitude up to 2.5e 2 /h call for an alternative ex- 
planation. In disordered graphene, close to the Dirac point, the mean carrier density is low 
and screening is reduced so that the surrounding disordered medium drives the system into 

formation of co-existing local 



a fluctuating potential landscape. This leads to the natura 
electron and hole puddles [l^]. Following the lines of 20, 



26], disorder can be generally 
described through a charged impurity density rii mp that scales 10 15 //i. Using 

mobility /i = 1300cm 2 jV.s for our sample, we compute n imp = 3.85 x 10 12 cm -2 , so that the 
typical mean scattering length is roughly 5nm. Using a scanning single-electron transistor, 
Martin et. al. [13] have demonstrated that the mean puddle dimensions are of the same 
order as the typical disorder length. It is therefore reasonable to consider a size distribution 
of electron-hole regions centered at £ p w 5 nm. When a magnetic field is applied to such 
an inhomogeneous system, within a specific puddle charge carriers tend to be quantized in 
cyclotron orbits provided the puddle size is larger than the magnetic length £b = J jj^- 
Therefore, as the magnetic field increases, the puddles successively enter into the quantum 
regime depending on their sizes. This process leads to conductance fluctuations of maxi- 
mum amplitude expected for £b = £ p i.e. B = 25T as shown in figures IH-a) and b). This 



27fl who analysed magneto- 



result should be compared with the work of Branchaud et. al. 
conductance fluctuations close to CNP and deduced a mean puddle size of about 125nm 
for clean graphene. Also, Rycerz et. al. 28| interpreted the conductance fluctuations as a 



natural consequence of the presence of electron-hole puddles, invoking a percolation regime 
in which the magnetic field induces small variations of the charge carrier trajectories. 

To conclude, we have investigated the electronic properties of disordered graphene in the 
close vicinity of the Dirac point. At low magnetic field, the system is well described using a 
simple two fluid model where both electron and hole-like particles contribute to transport. 
The co-existence energy interval is well accounted by existing theory for disordered graphene 
whereas the mobility dip at low carrier density is possibly explained by the presence of 
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vacancies. When the magnetic field is high enough so that the energy spectrum displays 
a clear LL quantization, large fluctuations of the magneto-resistance near the Dirac point 
are observed and explained in terms of topological transition of electron-hole puddles into 
the quantum regime. These fluctuations progressively vanishes as the magnetic field is 
further increased, suggesting a conduction regime where the magnetic field length is smaller 
than the mean puddle size. For clean graphene, spin and valley degeneracy splitting occur 
in strong magnetic field and, for the cleanest samples, the fractional quantum Hall effect 
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30| has even been observed. Many experimental studies converge to the assumption 
that the more disordered graphene is, the highest is the magnetic field at which such effects 
occur. It makes no doubt that similar experiments with samples having different degrees of 
disorder will significantly enlarge our understanding of the dirty regime. 

We are embedded to K.S. Novoselov for providing sample. This work was funded 
by the french National Agency for Research (ANR) through the program referenced 
ANR-08-JCJC-0034-01 and EUROMAGNET under E.U. contract n° 228043. 
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